Pseudogap in underdoped cuprates and spin-density-wave fluctuations 
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We analyze fermionic spectral lunction in the spin-density-wave (SDW) phase ol quasi-2D cuprates 
at small but finite T. We use a non-perturbative approach and sum up infinite series of thermal 
self-energy terms, keeping at each order nearly-divergent (T/J) \ log e| terms, where e is a deviation 
from a pure 2D, and neglecting regular T/J corrections. We show that, as SDW order decreases, the 
spectral function in the antinodal region acquires peak/hump structure: the coherent peak position 
scales with SDW order parameter, while the incoherent hump remains roughly at the same scale as 
at T = when SDW order is the strongest. We identify the hump with the pseudogap observed in 
ARPES and argue that the presence of coherent excitations at low energies gives rise to magneto- 
oscillations in an applied field. We show that the same peak/hump structure appears in the density 
of states and in optical conductivity. 

PACS numbers: 71.10.Hf, 75.10.Jm, 74.25.Dw 



I. INTRODUCTION 

Understanding of the phase diagram of cuprate super- 
conductors continue to be one of central topics in theo- 
retical condensed matter physics. 1 Parent compounds of 
cuprates are quasi-2D antiferromagnetic insulators, heav- 
ily overdoped cuprates are Fermi liquids. In between, 
systems are d— wave superconductors at low T < T c and 
display the pseudogap behavior at larger T c < T < T* . 
How an insulator transforms into a Fermi liquid and what 
is the origin of the pseudogap are still the subjects of in- 
tensive debates among researchers. 

The pseudogap region exists both in underdoped and 
overdoped cuprates, but the physics evolves substantially 
between these two limits. For overdoped cuprates, there 
is rather strong evidence^ that the pseudogap region is 
best described as a disordered superconductor, when the 
gap is already developed but the phase coherence is not 
yet set£~— In this doping range, fermions are reasonably 
well described as strongly interacting quasiparticles with 
a large, Luttinger-type underlying Fermi surface (FS). 7 
The d— wave pairing in this doping range most naturally 
originates from the exchange of overdamped collective 
bosonic excitations of which spin-fluctuation mediated 
pairing is the key candidate! 7 - - — In underdoped cuprates, 
situation is more complex. On one hand, ARPES data 
taken at low energies (below 50meV) and low T were in- 
terpreted as the indication that the underlying FS still 
has Luttinger form, and the gap extracted from the posi- 
tion of the still visible narrow peak in the spectral func- 
tion has a simple d— wave, cos 2<p form over the whole FS, 
including antinodal region around (0, it) and symmetry- 
related pointS i 12 i 13 On the other hand, ARPES data 
taken in the antinodal region show that the spectral func- 
tion in the pseudogap regime develops a broad maximum 
at around 100 - 200meVMr^ The jury is still outi£ 
whether the observed high-energy hump and low-energy 
peak are separate features, or the peak and the hump de- 
scribe the same gap, A(fc), which strongly deviates from 



cos 2(f) form with underdoping. The experimental results 
in Refs. [l^l - [l^ll7l - [20l were interpreted both ways. We 
side with the idea that the pairing gap remains cos 2<f> 
even in underdoped materials, and the hump is a separate 
feature, associated with Mott physics. We further take 
the point of view that the origin of the hump is the devel- 
opment of precursors to a Heisenberg-type antiferromag- 
netically ordered state at half-fillingi2 1 - ~ 24 ' 26 ~ — These pre- 
cursors are generally termed as SDW precursors though 
one should keep in mind that the half-filled state is the 
strong coupling version of SDW and is best described 
by the Heisenberg model with short-range exchange in- 
teraction. The SDW precursor scenario has been wildly 
discussed in mid-90th , 21 i 23 ~ 25 i 30 and is near-universally 
accepted scenario for electron-doped cuprates 24 ' 26 i 45 ' 46 . 
For hole-doped cuprates, it was, however, put aside for a 
number of years in favor of non- Fermi liquid type scenar- 
ios^ The SDW scenario, however, re-gained support in 
the last few years, after magneto-oscillation experiments 
in a field of 30 — 60T detected long-lived Fermi liquid 
quasiparticles near small electron and hole FSs<22 Such 
FS geometry is expected for an SDW ordered stated and 
early theory prediction was that a field drives the system 
towards a SDW instability^ 3 - Long-range antiferromag- 
netic order in applied field has been explicitly detected 
in recent neutron-scattering experiments on underdoped 
YBCO (Ref. H3). [Another widely discussed scenario of 
quantum oscillations, which we will not consider here, is 
a d— wave density- wave order.—] 

In this communication, we analyze the consistency be- 
tween the description of quantum oscillations and the 
pseudogap in underdoped cuprates within SDW scenario. 
The problem is the following: to explain quantum oscil- 
lations one has to assume the existence of small electron 
pocketsi 36 i 37 Such pockets do exist in the SDW scenario 
near (0,7r) and symmetrey-related points, but they are 
present only if SDW order (S Z (Q)} = (S z ) is smaller than 
a threshold (Q = (n, it)). For larger (S z ), only hole pock- 
ets around (tt/2,tt/2) are present, while excitations near 



2 



(0, 7r) have a gap of order At' ~ 0.2eP (see Fig. [2]). Antin- 
odal pseudogap detected in ARPES experiments in zero 
field is of the same magnitude. 13 A held of 40 — 60T is too 
small to affect energies of 0.2eV, hence the same 200meV 
pseudogap should be present in the ordered SDW stated. 
If this pseudogap is viewed as a precursor to SDW, one 
could expect that it simply sharpens up in the ordered 
SDW phase and transforms into the true antinodal gap. 
But then there will be no electron pockets in the SDW 
phase, in disagreement with magneto-oscillation experi- 
ments. The pseudogap and quantum oscillations can be 
reconciled within SDW scenario only if the evolution of 
the spectral function between paramagnetic and SDW 
states is more complex than just the sharpening of the 
pseudogap, and antinodal spectral function in the SDW 
phase contains both, high-energy pseudogap and electron 
pockets. This co-existence also explains ARPES data in 
the superconducting stat o 12 d 3 because pairing of coher- 
ent fermions near electron pockets gives rise to a sharp 
peak in the spectral function at the gap energy. 

To address this issue, we consider how fermionic 
Green's function G(k, ui) evolves within SDW phase, as 
SDW order gets smaller. We depart from Heisenberg 
antiferromagnet with exchange interaction J = At 2 /U 
and consider analytically how G(k,uj) is affected by ther- 
mal fluctuations which in quasi-2D systems destroy long- 
range order already at T <C J. We neglect regular T j J 
corrections but sum up infinite series of self-energy terms 
which contain powers of j3 = (T/ttJ)\ log e| , where e is 
a parameter which measures deviations from pure two- 
dimensionality and which we will just use as a lower cutoff 
of logarithmically divergent 2D integrals. In practice, e 
is e.g., the ratio of the hoppings along z-axis and in xy 
plane 3 ^. SDW order disappears when (3 = (3 cr = 0(1). 
This yields a set of integral equations for the Green's 
function and (S z ) which we obtain and solve. 

In the terminology of Ref. ScJ our computations are 
valid in the renormalized-classical regime of quasi-2D sys- 
tems. The idea that, in this regime, at low enough T, 
one cannot restrict with Eliashberg or FLEX approxima- 
tions and has to include self-energy and vertex correc- 
tions on equal footings has been put forward in Refsi 21 ' 23 . 
The computational procedure that we are using is sim- 
ilar to eikonal approximation in the scattering theory. 
Such procedure has been used in the study ID CDW 
systems by Sadovskii^i and others^ 2 -, and has been ap- 
plied to cuprates in Ref. [H to analyze SDW precursors 
in the paramagnetic phase (for latest developments, see 
Ref. 1271 ). Our computation has one advantage over ear- 
lier works: in the SDW-ordered state we don't need to 
assume that T is larger than some threshold To to re- 
strict with only thermal fluctuations (i.e., with the con- 
tributions from zero Matsubara frequency). All we need 
is a small e such that T/J\ loge| = 0(1) even when T/J 
is small. In a paramagnetic phase, eikonal approxima- 
tion is only valid when T > To, and To increases as one 
moves away from the SDW phase. We assume that near 
SDW boundary Tq is small and apply our theory also to 
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FIG. 1: (Color online) Schematic phase diagram of hole and 
electron-doped cuprates. The regions of antiferromagnetism 
(AF), superconductivity (SC) and pairing fluctuations (PF) 
are shaded. In the PF region, there are vortex excitations 
and large Nernst signal.— Precursors to SDW appear at a 
non-zero T due to strong thermal fluctuations. In the region 
where SDW precursors are already developed, the onset tem- 
perature for the pairing increases with increasing x (Ref. IFlh. 
in the region with no SDW precursors, it decreases with in- 
creasing x (Ref.0), the maximum T c is in the area between 
the two regimes. A similar phase diagram has been proposed 
in Ref. 

a paramagnetic phase. Our results for a paramagnet are 
in full agreement with Ref. [22|. 

Note that in our theory (and in Ref. l22t ). the para- 
magnetic state is a Fermi liquid at the lowest energies 
at T = 0. We don't discuss here a possibility that 
a new, non-Fermi liquid state emerges near the region 
where SDW order is lost^ 3 -. We also don't discuss possi- 
bilities of more complex spin order and of open electron 
Fermi surface o 36 i . Our results are expected to survive if 
SDW order is incommensurate (with Q still near (7r,7r), 
but whether our results survive if the system develops a 
stripe order remains to be seen. 

We found that the spectral function A(k,uj) = 
(l/ir)\ImG{k, uj)\ near (0, tt) in the SDW state has a peak 
and a hump. Both originate from a single peak at the 
value of the T = SDW gap at (0, tt). The hump moves 
little as SDW order decreases and just gets broader, while 
the peak follows (S z ), shifts to lower energies as SDW 
order decreases, and vanishes /? = /3 cr , when the sys- 
tem enter the paramagnetic phase. At f3 > f3 cr , only the 
hump (the pseudogap) remains, and the spectral func- 
tion at antinodal k = hp has camel- like structure, with a 
minimum at uj = 0. As /3 increases further, A(kF, uj = 0) 
increases and eventually the spectral function at k = hp 
develops a single-peak at uj = 0, as it should be for a 
system with a large, Luttinger FS. 

Rewinding this backwards, from a paramagnet to an 
SDW state, we see that the system first develops a pseu- 
dogap as a precursor to SDW. When SDW order sets in, 
the pseudogap sharpens up, but, in addition, there also 
appears a true quasiparticle peak at low-energies. The 
residue of the peak increases as (S z ) increases. When 
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(S z ) is below the threshold, electron pockets are present, 
and the spectral function near (0, w) has a low-energy co- 
herent peak and a hump at about the same energy as the 
pseudogap in a paramagnetic phase. When SDW order 
gets larger, electron pockets eventually disappear, peak 
and hump come closer to each other and merge when 
(S z ) reaches its maximum. 

This peak/hump structure also shows up in the den- 
sity of states and in the optical conductivity <j(uS). In 
the Mott-Heisenberg limit (2U(S Z ) is larger than free- 
fermion bandwidth) the conductivity at T — is zero 
up to a charge-transfer gap U ~ 1.7eV. Once SDW or- 
der gets smaller, the peak at U splits into a hump which 
slowly shifts to a higher frequency, and a peak whose en- 
ergy scales as 2U(S Z ). In addition, there appears a metal- 
lic Drude component at the smallest frequencies. This 
behavior is quite consistent with the measured <j{uj) in 
electron-doped cuprates, where SDW phase extends over 
a substantial doping range.— ~— 

We also found that, at a finite T, the system in 
the pseudogap phase retains the memory about pock- 
ets. There are no real pockets in the sense that there 
is no two-peak structure of the spectral function at zero 
frequency along zone diagonal, but we found that, when 
P > Per, the spectral weight at u> = extends almost all 
the way between the original FS at Uf and the "shadow" 
FS at k — Q — kp (see Fig. [5J. As f3 becomes larger, 
the k— range where the spectral weight is finite shrinks, 
and at large /3 the spectral function recovers Drude-like 
structure typical for a metal with a large, Luttinger FS. 

This analysis can be extended into a superconduct- 
ing state. A system does not need to possess coherent 
quasiparticles to develop a pairing instability 48 ! 49 but 
fcrmionic coherence emerges below the actual T c much 
in the same way as it emerges in the SDW ordered state. 
The spectral function in the antinodal region then dis- 
plays a coherent superconducting peak and a hump cen- 
tered at, roughly, the energy of the antinodal SDW gap 
at T = 0. This picture is consistent with the data from 
Refs. [Mil 

The overall conclusion of our analysis is the phase di- 
agram of the cuprates presented in Fig. fT) A similar 
phase diagram has been proposed in Ref. l5Cl At T ^ 0, 
there is a region where the system displays SDW pre- 
cursors. In this region, magnetic excitations are moder- 
ately damped, propagating magnons, and magnetically- 
mediated d— wave pairing interaction decreases as doping 
decreases due to a reduction of the electron-magnon ver- 
texj£i~— In this region, the antinodal pseudogap, caused 
by SDW precursors, and a d— wave pairing gap co-exist. 
Outside this region, SDW precursors do not emerge, and 
arcs and other pseudogap features are caused by thermal 
fluctuations of a pairing gap^ ' 17 ' 55 

We discuss the computational procedure in the next 
section, present the results in Sec. 3 and conclusions in 
Sec. 4. 
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FIG. 2: (Color online) Evolution of the FS with increasing 
SDW order, (a) - paramagnetic phase; (b) - SDW order is 
about to develop ({S z ) = + ). The shadow FS emerges, but 
the residue of fermionic excitations at the shadow FS is + ; (c) 

- small (S z ), both hole and electron pockets are present; (d) 

- a larger (S z ), only hole pockets around (7r/2,7r/2) remain. 



II. COMPUTATIONAL PROCEDURE 



Our point of departure is the mean-field, SDW the- 
or y 2 1 ' 57 ~ — of the antiferromagnetically ordered state in 
the large U quasi-2D Hubbard model at T — 0. We as- 
sume t — t' dispersion in the XY plane [e k — — 2t (cos k x + 
cos k y ) — At' cos k x cos k y ] and weak dispersion along the 
Z axis, which we will not keep explicitly in the formu- 
las. Mean-field description neglects quantum fluctuations 
and is rigorously justified when the model is extended to 
25 3> 1 fcrmionic flavors^ but qualitatively it remains 
valid even for S — 1/2 ((S z ) becomes 0.32 instead of 0.5). 

Long-range antiferromagnetic order splits the 
fermionic dispersion into valence and conduction 
bands, separated by U, and gives rise to a two-pole 
structure of the bare fermionic Green's function: 



Go (Lo,k)=u 2 k G c + v 2 k G^ 



(1) 



where u kl v k = T Atj k /E k )/2, j k 

cosfc y )/2, and 



(cos k x 



Gn — 



— 



UJ ■ 



El 



(2) 



The dispersions of conduction and valence electrons are 
given by E% v = ±E k 



At' cos k x cos k y — [i, where E k 



16* 2 7fc] 1/2 , A = U(S Z ), and fi « -A is the 



chemical potential. 

The shape of the FS depends on the value of Ao. We 
show the evolution of the FS with increasing A in Fig. [21 
For small Aq both hole and electron pockets are present, 
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for larger Ao only hole pockets remain. At large U/t, 
which we assume to hold, valence and conduction bands 
are well separated near half- filling at T = 0, u 2 v k « 
1/2, and the FS only contains hole pockets. 

The value of (S z ) is determined by the self-consistency 
condition 



(S.) = 



d 2 k 
(2-k)- 



:UkVk 



J ^n F (cj)Im [G c - Gg] , 



(3) 



where both G§ and Gq are retarded functions and uf(uj) 
is the Fermi function. At large U, and near half-filling, 
(S z ) ~ 1/2, and A [7/2. For larger dopings and 
smaller U, (S z ) is smaller already at T = 0. 

Fermion-fermion interactions in the ordered SDW state 
can be cast into interactions between fermions and 
magnons. These interactions are described by the ef- 
fective Hamiltonian^i 



el— mag 



= J2J2[ a ak a Pk+ae+V aa (k,q) 

oc,/3 k,q 

+ b lkhk+ q eqV hb (k,q) + a+ k b pk+q e+V ab {k,q) 



+ b lk a Pk+ q eqV ba (k,q) +h.c.]5 a -p, 



(4) 



where operators a and b represent conduction and valence 
band fermions, respectively. The vertex functions are 
given by 

V aa ,bb(k, q) = —= [±(u k u k+q - V k Vk+q)r] q 

+ (u k V k + q - V k Uk+q)fi q , 

V a b,ba(k,q) = -^=[(u k v k+q + v k u k + q )r)q 

T(u k u k+q + v k v k+q )fj q ], (5) 



with 



G = 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ • • • 
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FIG. 4: Diagrammatic series for the Green's function. Only 
thermal contributions in which a fermion jumps from a va- 
lence to a conduction band (and vise versa) and emits (ab- 
sorbes) a transverse spin-wave are included. 



A. One-loop perturbation theory 

There are several contributions to fermionic self-energy 
E(fc,w) to one loop order, but the earlier study by Morr 
and one of us has found^i that the dominant one at the 
lowest T comes from the interaction between valence and 
conduction fermions mediated by the exchange of low- 
energy transverse spin- waves (see Ref. [2l| for details) . To 
one-loop order, spin-wave mediated interaction gives rise 
to 



£?' c (w,k) =/3A 2 Q G c , v (w,k) 



P_3 



(7) 



where, we remind, /3 — (T/irJ)\ loge|. The order pa- 
rameter also acquires a correction proportional to j3. To 
obtain it, one has to substitute the self-energy into the 
Green's function and compute (S z ) using Eq. ([3]), but 
with the full G instead of Gq. This yields 



(8) 



Vq 



1 /1-79X 1 / 4 



V2A1 + 7, 



1 / 1+7^1/4 f 



At the mean-field level, the spectral function A(k, ui) 
has two S— functional peaks at u = Ef. and uj = E^, and 
the density of states has a gap 2Ao = U. Our goal is to 
analyze how this spectral function gets modified once we 
use Eq. (|4]) and compute thermal fermionic self-energy 
and thermal corrections to (S z ). 




(a) 



(b) 



(c) 



FIG. 3: Three equivalent diagrams for two-loop corrections 
to the Green's function. The second and the third diagrams 
contribute to fermionic self-energy. 



i.e., A 2 = Aq(1 - /3 + 0((3 2 )). The 0((3) (| log e|) correc- 
tion to (S z ) is in agreement with Mermin- Wagner the- 
orem. However, when we combine self-energies (J7J and 
Gq' c in which A = U/2 is replaced by A = U(S Z ) and 
obtain the new Green's function, we find that 0(/3) terms 
cancel out, i.e., to first order in j3 the fermionic Green's 
function does not change: 



G(w,k,0) 



1 



1 



u - E c - SJ 



- 16t 2 7 2 - (A 2 + /3A 2 ) 



oj - E v - 
= Go(w,fc), (9) 



where u> = lu + At' cos k x cos k y + ji. This result was ob- 
tained ir>2i and was interpreted as an indication that the 
SDW form of G(fc, uj) may survive even when (S z ) van- 
ishes. However, one-loop result is at best indicative, and 
we need to go to higher orders to verify what happens 
with the fermionic Green's function when ft increases. 
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B. Two- loop corrections 

As the next step, we obtain two-loop formulas for the 
self-energy and (S z ). The two-loop diagrams for the self- 
energy are the second and third diagrams in Fig. [3l Eval- 
uating them in the same approximation as one-loop dia- 
gram, we obtain 



(w,k) = 



2/3 2 A 4 c 



( W - m) (uj - E?y 



£2 (w,k) = 



(w - £^ (w - 



(10) 



Substituting these self-energies into the valence and con- 
duction Green's functions together with one-loop dia- 
grams, we obtain after a simple algebra 



G(w,k) 



16i 2 7 2 - (A 2 + 2A 2 /3) 



16t 2 7 2 - (A 2 



(11) 



Evaluating (S^) in the same two-loop approximation we 
obtain 



1 



5fP 
2 8 



(12) 



a.u 



i • 1 — 

A (co) 




A 7 (co) 




a.u. 



FIG. 5: (Color online) The spectral function A(u , khs , /3) at a 
hot spot. Here and below a.u. stand for arbitary units. Left 
panel - Atr, obtained by including only transverse spin waves. 
Right panel - Ai+tr, obtained by treating transverse and lon- 
gitudinal spin excitations on equal footings. Dashed lines - 
f3 — 0.6f3 cr , solid lines - ft = 0.8/3 cr . A tr has a branch cut at 
energy which scales with the magnitude of SDW order param- 
eter, while Ai+tr has a hump at energy which roughly remains 
the same as SDW gap at T = 0. The actual A(u,khs, @) at 
f3 < /3 cr coincides with Atr at low frequencies, and crosses 
over to Ai+tr at frequencies larger than the gap for longitudi- 
nal fluctuations. As a result, the actual spectral function has 
a peak at a low energy and a hump at a higher energy. 



such that A 2 = Ag (l - (3 + 3/3 2 /2 + 0{/3 3 )). Substitut- 
ing now this A 2 into (|TT|) we find that, up to two-loop 
order, 



GV,k) = - 



1 



uj- 



16t 2 7 2 



A? 



16i 2 7 2 - A 2 



(13) 



where Ai = A (l + 0/2 + 5/3 2 /8) and A 2 = A (l - + 
/? 2 /4). 

We see that the Green's function splits into two com- 
ponents. Both have SDW form, but the values of A are 
different - A2 decreases with while Ai increases. This 
implies that the peak in the spectral function, originally 
located at ui 2 = 16i 2 7 2 + A§, splits into two subpeaks - 
one shifts to higher |a)|, another to smaller |cD|. This is the 
new trend, not present in the one-loop approximation. 

This consideration also shows that to understand what 
happens when — O(l), one cannot restrict with a few 
first orders in the loop explansion but rather has to sum 
up infinite number of terms. This is what we are going 
to do next. 



• the SDW order parameter (S z ) is given by the same 
loop expression as in the mean-field theory, Eq. ^ , 
but with the full G instead of Go 



in turn appears in the 

V ,c 



the renormalized (S 

Green's function through Gq' c in which Ao 
has to be replaced by A = U(S, 



U/2 



• the full Green's function G(w,k, 0) is given by 
G(w,k,/3) = v 2 k G v (oj : k,P) +u 2 G c (u;,k,/3), where 
G v,c are the full Green's functions for valence and 
conduction fermions and Uk and v% are the same as 
in ([1]) but with A instead of Ao 

• at loop order n of the perturbation theory there are 
(2n — 1)!! equivalent contributions to the full G c,v , 
each contains (0G%(u, k)Gg(w, k))" (see Fig. Q). 

Because of the last item, it is advantageous to sum up 
infinite series of diagrams for the Green's function rather 
than for the self-energy. We have 



G?;>,k,/3) 



(w,k) + ^(2n- 



l)!!GS'°( W ,k) 



C. Non-perturbative Green's function 

We list several results which can be explicitly verified 
by doing loop expansion order by order in 0: 



/3A 2 G V,k)G§( W ,k) 



(14) 



where dots stand for non-logarithmic corrections and 
subindex tr implies that we only considered interaction 
with transverse spin waves. 
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Substituting the expressions for Gq' c and summing up 
asymptotic series we obtain 



1/2 



x exp 



^ 2 - g| 



A \2p) (Co 2 ~E 2 y/ 2 

{ r /^2~ 

<^ i + Erfi 



(15) 



I to 2 - El 
2A 2 p 



where Erfi(z) = — iErf(iz) is imaginary error function 
(Erfi(ir) is real when x is real and imaginary when x is 
imaginary). Observe that ImG\f vanishes when 
El 

To one- loop order, Eq. (fTS")) reduces to Gq' c , but be- 
yond one loop the Green's functions G v,c obviously be- 
come /3-dependent. The spectral function A tr (uj, k, 0) = 
ir~ 1 \ImG(oj — iO, k, /3)|, is readily obtained from fp~5|) : 



A tr (u,k,0) = 



A 

■«2 



i 



(c^ 2 - £2)1/2 



(16) 



where = 1 for x > and zero otherwise. 

Substituting G"' c from (jT5]) into the expression 
for (S z ) we obtain how the SDW order parameter 
evolves with 0. 



(S z 



d 2 k 
(2tt) 2 



2 A 



du>- 



exp 



x tlf(ljj — /i — 4t' cos k x cos ky)0 oj 



2 -El 
El], 



(17) 



The remaining unknown parameter /i is fixed by the con- 
dition on the number of particles in the SDW stated: 



(1-x) 



cPkdcu 

jA tr {u, k, p)riF{u). 



(2nf 



(18) 



These coupled equations were solved numerically. The 
dependence of (S z ) on /3 or doping x is quite as ex- 
pected: (S z ) monotonically decreases as /3 or x increase, 
and vanishes at some particular p cr and x cr (see Fig. [6]). 
The spectral function A tr (uJ, k, 0) has sharp 8— functional 
peaks at /3 -> 0, at Co = ±E k (uj = E v k '°). At fi- 
nite quasiparticle peaks transform into branch cuts at 
\Co\ = Ek, with the width of order /3 cx T, and the spectral 
weight extends to larger frequencies. Near the branch cut 
A tr (w,k, 0) diverges as l/y^aL We plot A tr (w,k,j8) at a 
hot spot k hs = (k x , 7T - fc x ) in Fig. [5] 



D. Further modifications of the spectral function 

On a more careful look, we found that the spectral 
function given by Eq. (I16|) have to be further modified 
by two reasons. First, in the calculations above we only 



included the self-energy due to exchange of transverse 
spin waves, and neglected the self-energy due to exchange 
of longitudinal spin fluctuations. This is justified at small 
when longitudinal fluctuations are gapped, but when 
P ~ /3 cr , longitudinal fluctuations are nearly gapless and 
are as important as transverse ones. The limiting case 
when transverse and longitudinal spin propagators are 
identical can be studied within the same approximation 
as before, the only difference is that combinatoric factors 
are now (2n + l)!!/2" (Ref. [H). As a result, the spectral 
function becomes 



A i+tr (u;,k, 0) 



ul\u) + E k \ 



1 




exp< 



-(a; 2 



El) 



El)- 1 ' 2 



Co 2 - E 2 



(19) 



where subindex I + tr implies that this is a contribution 
from both longitudinal and transverse spin excitations. 
One can easily make sure that Ai +tr is also fully inco- 
herent at P > 0, but, in distinction to A ir , it vanishes at 
Cu = ±E k and has a hump at a frequency which remains 
of order A for all j3 < (3 cr We plot A/ +ir (cj, k, 0) at a 
hot spot in Fig. El The function Ai +tr (ui, khs,0) becomes 
particularly simple at /3 > p~ cr : 



Aj +tr (w,khg) 



: exp-; 



0A§ 



(20) 



For a generic /3 < /3 cr , Eq. (fT9"f is the correct result 
at high energies, larger than the longitudinal gap, and 
Eq. (p~6|) is the correct result at smaller energies. As /3 
approaches p~ cr the range of applicability of A tr (u), k/j S ) 
shrinks. The full formula cannot be obtained within 
eikonal approximation, but it is clear that the actual 
A(w,k, 0) contains both, the branch cut, l/y/x singu- 
larity near uj = E k ' c and the hump at a frequency of 
order Ao, where the quasiparticle peaks were located 
at /3 = 0. To simplify the computational procedure, in 
Fig. IH1 below we use for the actual A(oj, kh s ,0) the func- 
tion A tr (ui, khs, 0) up to a frequency where it crosses with 
Ai +t r, and use the function A; +tr (u;, khs, 0) at larger fre- 
quencies. 

Second, Eqs. (fT()]) and (fH)j) show that fermionic coher- 
ence is lost immediately when /3 becomes non-zero (the 
pole transforms into a branch cut). Meanwhile, from 
physics perspective, as long as the system has a SDW or- 
der, a Fermi liquid behavior near the pocketed FS should 
be preserved, i.e., a quasiparticle peak with T 2 logT 
width should survive, albeit with a reduced magnitude. 
The reason it was lost in the calculations above is because 
we completely neglected regular classical and quantum 
corrections to the Green's function (dots in Eq. ([13])). 
We verified that, when these terms are included, only a 
part of Gq'° gets involved in the renormalizations by se- 
ries of /3 n corrections, the other stays intact. This implies 
that the the actual A v £ n = Z p A" a :c + (1 - Zp)A v < c . The 
residue Zp is some number 0<Z J g<lat/3 = 0, where 
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FIG. 7: (Color online) The appearance and evolution of the electron pocket near (k x , k y ) = (0, n) and symmetry related points. 
Electron pocket appears as a single point at (0, n) once /3 increases and reaches a critical value, and evolves with increasing /3. 
From left to right: p//3 cr = 0.8; 0.85; 0.95. We set x = 0.02. 
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anyway A v > c — A V q C , it decreases as /3 increases, and van- 
ishes at /3 — (3 cr . For definiteness, we used Zp = 0.2 in 
the panel < (3 < f3 cr in Figs. HJandlU 



FIG. 6: (Color online) The SDW order parameter (S z ) vs /3 
at a given x = 0.05. For simplicity, we set Zp to be a constant 
(= 0.2). Inset: (S z ) vs x at a given ft = 1. 



These two additions also affect the formula for (S z 
and the equation for which become 



J 



and 



I S^^[n F {El)-n F {Ei)] 



(2ny E k 



d 2 k 
W) 2 



,T"A f°° _ ex P 
— / du) 



f a 2 -El \ 



El 



up (cj — fj, — At' cos k x cos k y 



El 



(1 



du 



_d 2 k 

(2 



(21) 



2 -!Ln F (u,) x | ^ [t%6(w - El) + u\8{w - E c k )} + (1 - Z P )A( U , k)| . (22) 



III. THE RESULTS in FigO We used several phenomenological forms of Zp 

in which Z$ CT — 0, but found that the functional forms 
of the spectral functions do not depend on Zp in any 
We solved Eqs. ([21]) and ((22]) numerically and plot the substantial way. For simplicity, and because below we 
dependence of the order order parameter (S z ) on /3 and x only present the results for a few /3 in the SDW state, 
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FIG. 8: (Color online) The spectral function, Af u tl(w, k, 0) at ui = along the diagonal (nodal) direction in the Brillouin zone. 
We used Z — 0.2 for the top right panel. Observe that the systems retains a memory about a "shadow" FS when SDW order 
disappears at = /3 cr (the spectral weight is non-zero in the whole region between the original and the shadow FSs). We set 
x — 0.05, t = 0.32Aq, t' — —0.2t. For these parameters, /3 cr ~ 1. 



we set Zp — 0.2, independent on the actual (3. Also, all 
cases (even when Zp = 0) (S z ) monotonically decreases 
when either x or f3 increase, and vanishes at some x cr 
and (3 cr . 

In Fig. [7] we show the evolution of the electron Fermi 
surface near (0, ir) and symmetry-related points. The 
electron pocket is absent at small j3, when SDW order is 
strong, but appears when (3 exceeds some critical value. 
The electron pocket evolves with increasing f3 and even- 
tually disappers when the system looses long-range SDW 
order. This behavior is qualitatively consistent with the 
mean-field SDW picture. 

The results for Af u u(oj, fc,/3) are presented in Figs. El 
HI andfTUl In Fig. [5] we plot the full spectral function at 
zero frequency along the diagonal direction in the Bril- 
louin zone. In the two limits, (3 = and (3 3> /3 cr , 
the system possesses sharp quasiparticles - in the first 
case at the two sides of the hole pocket, in the second 
case at the large, Luttinger FS. In between, the spectral 
function evolves, as j3 increases, from a well-pronounced 
two-peak structure to a completely incoherent structure 
at j3 = /3 cr , in which the spectral weight at ui = is 
spreaded between the original and the "shadow" FSs (in 
reality, Af u u(0, fc, j3) spreads outside of this range, but 



to find these tails of the spectral function one has to go 
beyond the accuracy of our calculations) . This result im- 
plies that the system does retain some memory about 
SDW pockets even when (3 = (3 cr and A = 0. When /3 
becomes larger than f3 cr and SDW order disappears, the 
region where Af u u(0, k, /3) ^ progressively shrinks with 
increasing f3 towards a single quasiparticle peak. We em- 
phasize that this evolution of Af u u{0, k, (3) with (3 is very 
different from that in the mean-field SDW theory, where 
the spectral function in the SDW state has two peaks at 
the two sides of the hole pocket, and the "shadow" peak 
just disappears when A vanishes. 

In Fig. [9] we show the frequency dependence of 
Af u u(u},khs, (3) at a hot spot for a wide range of f3. In 
the two limits, the behavior is again coherent: there are 
quasiparticle peaks at ui = w + — 4i'cos 2 k x = ±Ao deep 
in the SDW phase, and the peak centered at uj = oj = 
deep in the normal phase. In between, the spectral func- 
tion is again predominantly incoherent. Specifically, as f3 
increases towards (3 cr , the quasiparticle peak splits into 
the peak at Co — A, and the hump at a larger frequency. 
The frequency where the peak is located shifts down- 
wards with decreasing /3, while the hump remains roughly 
at around Aq, and, on a more careful look, shifts towards 
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FIG. 9: (Color online) The spectral function Af u u(co, kf, s , /3) at a hot spot, at various j3. The frequency is in units tj/Ao, where 
uj — u) + fi — At' cos 2 fc^.The parameters are the same as in Fig. [8] 



somewhat larger frequency. We show this behavior in 
more detail in FigfTUl where we plot the spectral func- 
tion for a range of (3 < f3 cr . The peak in the spectral 
function is the property of At r and Ao, and the hump 
is the property of Ai +tr - We emphasize that the peak 
and the hump do co-exist in the SDW phase, the first 
describes coherent low-energy excitations, the second de- 
scribes fully incoherent high-energy excitations. Once 
SDW order disappears at /3 cr , coherent excitations also 
disappear, but the hump remains and disappears only at 
much larger /3. It is quite natural to identify the hump 
at /3 > j3 cr with the pseudogap, and low-energy coherent 
excitations existing at f3 < f3 cr with the building blocks 
for magneto-oscillations. The implication of this result 
is that the 200meV pseudogap observed in ARPES in 
zero fielcUs is not an obstacle for observing magneto- 
oscillations once the system becomes SDW-ordered. 



A. Density of states 

The density of states (DOS) in the SDW phase, 
N(u>,P) = J(d 2 fc/47r 2 )A(w,k,/3), is obtained from 
Eqs. (JJ, (HU), and JTJ]). The expressions for N(u>) at 
finite P are rather complex for t' ^ but are simplified 
for t 1 — which we assume to hold in this subsection. A 



finite t' affects the behavior of the DOS at the smallest 
frequencies, but not at frequencies u> + yU > A, which we 
are chiefly interested in. 

For free fermions, we then have 



N (u,/3) 



dudv 



-l \/l — u 2 \J\ — v 2 



(23) 

At 2 (u + v) 2 ) 



2ir 2 ty/ti T 



A 2 



--K 



1 - 



A 2 



m 2 



where K(x) is the complete elliptic integral of the first 
kind (defined as in Ref. 62), and for t' = 0, u> = u> + /i. 
The DOS vanishes at |w| < A, diverges as l/\/x at \ui\ — 
A + 0, monotonically decreases at larger frequencies, and 
discontinuously drops to zero at the bandwidth, when 
\ui\ = VA 2 + 16i 2 . 

The DOS N tr (uj,f3) obtained using the spectral func- 
tion A tr from (1TB)) reduces to a ID integral 



JV tr (w,/3) 



2tt 5 /2^ 



dze 



xK 



A 2 - 2(3A 2 z 2 
16t 2 ~ 



2f3A 2 z 2 
(24) 



10 




a.u 




4 -2 2 4 



ml A„ 



FIG. 10: (Color online) The spectral function Af u ii(u), kh s , ft) 
at a hot spot, at various < f3 < f3 cr . The lowest plot 
corresponds to j3 — f3 cr , values of /3 decrease from bottom up. 
The frequency is in units uj / Ao . The parameters are the same 
as in Fig. [H] As /3 decreases, the peak and the hump come 
closer to each other. 



for < uj 2 - A 2 < 16t 2 , and 



N tr (u,j3) 



2n 5 / 2 t 



dze 



^JCj2 _ A 2 _ 2 (3A 2 z 2 



xK 



A 2 - 2 (3A 2 z 2 
16i 2 



(25) 



for <I> 2 — A 2 > 16i 2 . [The trick how to integrate over 
2D momenta k x and k y in J d 2 k/(4TT 2 )A tr (uj, k, /3) is to 
introduce u = cos k x and v 
the identity 



cos k y as new variables, use 



El 



2J2/3A 2 



dze~ z 8(u> 2 -E 2 k -2pA 2 z 2 ), 
(26) 

and use 5— function to perform the momentum integra- 
tion]. The density of states N tr vanishes at \uj\ < A, di- 
verges logarithmically at \uj\ = A + 0, and monotonically 
decreases at larger frequencies. There is no threshold 
at the bandwidth as At r (o>,k, /3) is nonzero everywhere 
at | a) | > A, but, indeed, at frequencies larger than the 
bandwidth our approximation eventually breaks down. 

For the spectral function Ai +tr from (fTT?)) , we obtain, 
using the same trick, 



ir 5 / 2 t 



xK 



A 2 - ^A 2 z 2 
16t 2 



f3A 2 z 2 
(27) 



FIG. 11: Density of states vs cj/Ao for different values of f3. 
We set t' = and t = 0.32A . 



for < ui 2 - A 2 < 16i 2 , and 



0J 



dzz e 



;,2 _ a2 _-, f> t 2 



•£ Vw 2 - A 2 - ^A 2 z 2 



xK 



_ uj 2 - A 2 - (3A 2 z 2 

1 m 2 



(28) 



for uj 2 — A 2 > lQt 2 This spectral function is continuous 
at \ui\ — A and has a broad maximum at frequencies 
comparable to A . 

This behavior is very similar to the one for the spec- 
tral function. Again, the transverse-only contribution 
Ntr nas a peak at the frequency ui = A, which scales 
with the order parameter of the SDW phase, while Ni+tr 
which treats contributions from transverse and longitu- 
dinal fluctuations on equal footings (as if longitudinal 
excitations were massless) has no features at A, but has 
a broad maximum at a frequency which remains compa- 
rable to Ao. Just like the spectral function, the actual 
DOS N(u>, f3) interpolates between these two terms. At 
frequencies smaller than the gap for longitudinal spin ex- 
citations, N{uj,(3) fa N tr {uj,(3), while at larger frequen- 
cies N(oj,(3) gradually approaches Ni +tr (uj, (3). 

We plot the DOS in FigHU We again set N(uj,/3) = 
N tr (uj,f3) at frequencies smaller than the crossing point 
between N tr {uj,f3) and Ni+t r (uj, f3), and set N(uj,f3) = 
Ni +tr {uj, P) at higher frequencies. And we again assumed 
that the full DOS Nf u u(uj, f3) is the sum of the incoherent 
N(ui,(3) with the factor 1 — Zp and the coherent, mean- 
field No(ui,f3) with the factor Zp. 

The plots clearly show the same trends in the DOS 
as we just discussed. At small /3, there is a sharp gap 
2Ao « U between valence and conduction bands. As (3 
increases, sharp gap decreases, the spectral weight ex- 
tends to higher frequencies, and the DOS develops a 
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hump at an energy comparable to Ao (i.e., the distance 
between the humps remains U). The peak in the DOS 
at the boundary of the sharp gap gets smaller as j3 in- 
creases and the sharp gap decreases. At /? = j3 cr the gap 
and the peak disappear, and the spectral function only 
possesses a hump. As j3 increases even further, the DOS 
at ui — increases (and also \i gets reduced such that uj 
comes closer to uj), and eventually the DOS at low fre- 
quencies recovers weakly-frequency-dependent form of a 
metal with a large FS. 

This physics is somewhat spoiled in the plots of 
Nf u u(u), (3) by the change of the behavior of N tr +i(uJ, f3) 
at the bandwidth (at uj = \/A 2 + 16t 2 ). This change of 
behavior is seen in Fig. [IT] as the discontinuity in the fre- 
quency derivative of Nf u u(uj,f3). For large /? the hump 
in the DOS is predominantly the effect of the bandwidth. 
However, at at intermediate ft, the hump is located at a 
frequency below the bandwidth, as is clearly visible in the 
top right panel of Fig. [TTJ (f3 = 0.1/3 cr ), and is therefore 
due to the physics that we described above. 



B. Optical conductivity 

The peak/hump structure also shows up in the optical 
conductivity. We computed the conductivity by standard 
means: by convoluting two full spectral functions Af u u 
using Kubo formula: 



a.u 



a(oj) 



d 2 k 



np(Q) — np(fl + uj) 



A fuU (n, k, /3)A fun (fl + u, k, /3), 



UJ 



(29) 



where Vk is fermionic velocity. We used the same compu- 
tational procedure as before: combined mean-field con- 
tribution to A, with the prefactor Zp, and incoherent 
part of A with the prefactor 1 — Zp, and used A tr for 
the incoherent part at small frequencies and Ai +tr at 
high frequencies. The results are shown in Fig. [T^] We 
again set t' = to simplify the calculations. At small /3, 
the conductivity almost vanishes up to uj = 2Aq = U. 
When /3 gets larger, the discontinuity at U splits into a 
hump which slowly shifts to a larger energy and a peak 
which scales with A and shifts to a smaller frequency as 
SDW order gets weaker. In addition, there also appears 
a Drude component at the smallest frequencies. 

The hump at large f3 > f3 cr is predominantly the effect 
of the bandwidth. It appears quite "sharp" in the last 
panel in Fig. I12[ but this is the consequence of our ap- 
proximation in which we neglected all regular self-energy 
terms. When these terms are included, the hump should 
definitely get broader. 

This behavior of <j(uj) is quite consistent with the ob- 
served evolution of cr{uj) in the SDW phase of electron- 
doped Nd2~ x Ce x Cu04, where the region of SDW- 
ordered phase extends over a substantial doping range, 
up to x ~ 0.15. Onose et al^i and others^ have found 
that conductivity does have a peak/hump structure at 




GT/A, 



FIG. 12: Optical conductivity, cr(w), for different j3 in the 
SDW phase at x = 0.02. For simplicity, we used Zp — 0.12 in 
all plots. At P = cr(u>) has a true gap and discontinuity at 
the onset of scattering between conduction and valence bands. 
The onset frequency is at ui slightly smaller than 2Ao = U 
because x is non-zero and the actual gap A < Ao. At finite /3 
the discontinuity splits into a hump which slowly moves to a 
higher frequency and a peak which moves to a smaller energy. 
We verified that the peak position scales with {S z ). 



finite x, the hump shifts towards a higher frequency with 
increasing x (from a charge-transfer gap of 1.7 'eV at x = 
to over 2eV at x — 0.1), while the peak shifts downwards 
as x increases, and was argued^* to scale with the Neel 
temperature T^. These two features are reproduced in 
our theory. 



IV. CONCLUSIONS 

To summarize, in this paper we obtained fermionic 
spectral function, the density of states, and optical con- 
ductivity in the SDW phase of the cuprates at small but 
finite T. We adopted non-perturbative approach and 
summed up infinite series of thermal self-energy terms, 
keeping at each order nearly-divergent T/J\ loge| terms, 
where e is a deviation from a pure 2D, and neglecting 
regular T/J corrections. We found that, as SDW or- 
der decreases, the spectral function in the antinodal re- 
gion acquires peak/hump structure: the peak position 
scales with the SDW order parameter, while the inco- 
herent hump remains roughly at the same scale as at 
T = 0, when SDW order is the strongest. We identified 
the hump with the pseudogap observed in ARPES exper- 
iments and identified coherent, Fermi liquid excitations 
at low energies as building blocks for magneto-oscillations 
in an applied field. The same peak/hump structure ap- 
pears in the DOS and in the optical conductivity. The 
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gap in the DOS scales with the SDW order parameter and 
disappears when SDW order vanishes, however the DOS 
also develops a hump at an energy which remains close to 
U/2, no matter what is value of the SDW order parame- 
ter. Optical conductivity at finite j3 has a Drude peak at 
the lowest frequencies, a peak, which again scales with 
the SDW order parameter and moves to smaller frequen- 
cies as P increases, and a hump which remains roughly 
at u + fj, ~ U and slightly shifts to higher frequencies as 
SDW order gets weaker. 

A more generic result of our study is the phase diagram 
for the cuprates shown in Fig.[TJ A similar phase diagram 
has been proposed in Ref. [5fJ. At large enough hole dop- 
ing and small enough temperatures, thermal fluctuations 
are weak and no SDW precursors appear. In this region, 
FS is large, and the physics is governed by the inter- 
action between fermions and Landau-overdamped spin 
fluctuations. This interaction gives rise to a fermionic 
self-energy which is Fermi liquid-like at the lowest ener- 
gies, has a non-Fermi liquid form, (iui) a , a < 1 at high 
energies, and displays a marginal Fermi liquid behavior in 
the cross-over region-. The same interaction with over- 
damped spin fluctuations gives rise to a d— wave pairing 
instability^*^ The onset temperature for the d— wave 
pairing increases as x decreases and approaches the uni- 
versal scale of around T p ~ 0.02vp/a, where vf ~ leV*a 
is the Fermi energy (Ref. |49( ) . The actual superconduct- 
ing T c is lower due to fluctuations of the pairing gap. 
In this regime, the thermal evolution of the FS is en- 
tirely due to thermal effects associated s— wave fermionic 
damping induced by scattering on thermal bosons) 5 -** 7 . In 
particular, at a finite T, the spectral function is peaked 
at zero frequency in some range of k around the nodal 
direction, despite that the pairing gap itself has cos 20 
form. 

At larger T and smaller x > 0, thermal fluctuations 
get stronger and give rise to SDW precursors. These 
does not imply that the FS actually becomes pocket-like, 
but the spectral function in the antinodal region develops 
a hump at a finite frequency, and the low energy spectral 
weight progressively fills in the area between the actual 
FS at k F and the "shadow" FS at k = (tt, vr) - k F . This 
behavior is at least qualitatively consistent with recent 
ARPES experiments aimed to verify the existence of the 



outer side of a pocket^ We also emphasize that, as long 
as SDW order is weak, the incoherent pseudogap peak 
is the dominant feature in the ARPES spectrum. This 
last observation is particularly relevant to LBCO near 
1/8 doping, for which there are indications of SDW or- 
dering. This SDW order is in any case smaller than the 
order at half-filling, and it is very likely that the domi- 
nant ARPES intensity is remains the preudogap despite 
the potential apearance of SDW order. This would be 
consistent with the observations in Refs. [20ll63l 

The redistribution of the fermionic spectral weight in 
turn affects the pairing problem: in the presence of SDW 
precursors, T c decreases with decreasing doping— be- 
cause now the pairing is due to the exchange of propagat- 
ing magnons rather than overdamped spin fluctuations, 
and the electron-magnon vertex gets smaller as SDW pre- 
cursors get stronger . 21 ' 52 i 53 A closely related effect which 
also reduces T c is the removal of the low-energy spectral 
weight due to pseudogap opening^ 4 - 

The precursors to SDW also give rise to dome-like 
behavior of the onset temperature for the pairing in 
electron-doped cuprates (left-hand side of Fig. [T]) but 
there the effect is weaker simply because pairing correla- 
tions are weaker! 47 ' 65 

Recently we became aware of a study of thermal SDW 
fluctuations by M. Khodas and A.M. Tsvelik (Ref^ 6 -). 
Their and our results are similar (e.g., both give linear 
in T width of the quasiparticle peaks at low T), but 
not identical as we studied isotropic quasi-2D syatems 
with a true long-range SDW order below a certain small 
T (P < 1) and exponential behavior of the correlation 
length at larger T, while Khodas and Tsvelik considered 
a 2D system with an easy-plane anisotropy and put spe- 
cial emphasis to the fact that spin correlations decay by 
a power-law at small T. 
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